arXiv: 1507.03140v3 [cond-mat.mes-hall] 11 Nov 2015 


Helicons in Weyl semimetals 


Francesco M.D. Pellegrino/’* Mikhail I. Katsnelson/ and Marco Polini^ 

^NEST, Scuola Normale Superiore, 1-56126 Pisa, Italy 
^Radboud University, Institute for Molecules and Materials, NL-6525 AJ Nijmegen, The Netherlands 
^Istituto Italiano di Tecnologia, Graphene Labs, Via Morego 30, T16163 Genova, Italy 

Helicons are transverse electromagnetic waves propagating in three-dimensional (3D) electron 
systems subject to a static magnetic field. In this work we present a theory of helicons propagating 
through a 3D Weyl semimetal. Our approach relies on the evaluation of the optical conductivity 
tensor from semiclassical Boltzmann transport theory, with the inclusion of certain Berry curvature 
corrections that have been neglected in the earlier literature (such as the one due to the orbital 
magnetic moment). We demonstrate that the axion term characterizing the electromagnetic response 
of Weyl semimetals dramatically alters the helicon dispersion with respect to that in non-topological 
metals. We also discuss axion-related anomalies that appear in the plasmon dispersion relation. 


Introduction. —After lightning strikes, very low- 
frequency transverse electromagnetic waves called 
“whistlers”^ propagate in the ionosphere, from one 
hemisphere to the other, along the Earth’s magnetic 
field lines. Interestingly, whistlers have a solid-state 
analog, which is usually called “helicons” . These 
transverse electromagnetic waves propagate in three- 
dimensional (3D) uncompensated metals subject to 
a uniform static magnetic field B = Bz. Helicons 
ultimately stem from the existence of the cyclotron 
resonance®, a single-particle excitation of the 3D 
electron system occuring at the cyclotron frequency 
uic = eB/(mc), where m (—e) is the electron mass 
(charge) and c is the speed of light in vacuum. In ordi¬ 
nary metals, helicons propagating along the magnetic 
field direction with a wave vector q = qz have the 
following free-particle-like dispersion relation®, 

^h(g) = 7 ^ < Wc , (1) 

zmh 

where mh = hujpj{2u!cC^) oc 1/B is the helicon effective 
mass and Wp = Airnee^/m the usual 3D plasmon fre¬ 
quency^’®. Helicons are relatively useless as a probe of 
many-body effects in the metallic medium but are inter¬ 
esting since they can hybridize with sound waves or dis¬ 
play interesting damping behaviors when the magnetic 
field direction is tilted from their propagation direction®. 

In this work we demonstrate that the helicon disper¬ 
sion in Weyl semimetals (WSMs)®^^® is greatly altered 
with respect to the textbook result (1). In particular, 
we show that helicon propagation in WSMs carries pre¬ 
cious information on the space- (b) and time-like (6o) 
components of the axion angle. We also highlight axion 
anomalies in the plasmon sector, which, to the best of 
our knowledge, have not yet been noticed. 

Maxwell equations in WSMs. —WSMs are recently dis- 
covered^^^^® 3D topological metals displaying an in¬ 
triguing electromagnetic response and Fermi-arc surface 
states. For the sake of simplicity, we here consider the 
Hamiltonian of a WSM with two nodes only^®: 

% = hvDT^cT ■ (—iV -f T^b) -b hr^bo . (2) 


Here, ud is the Dirac-Weyl velocity, describes the node 
degree of freedom with chirality ±1, and the 3D vector of 
Pauli matrices cr = (cr“, cr^, tr^)"'" describes conduction- 
and valence-band degrees of freedom. The two Weyl 
nodes are located at ±6 and shifted by 2 hbo in energy. 
Our results below can be easily generalized to the case of 
more than two Weyl nodes. Furthermore, in the present 
work we study helicons (and plasmons) in bulk Weyl 
semimetals, while the interplay between finite-thickness 
effects in a slab geometry (such as Fermi arcs) and elec¬ 
tromagnetic wave propagation will be the scope of future 
works. 

It has been demonstrated®® that the terms propor¬ 
tional to bo and b = {b^, by, bz)'^ in Eq. (2) can be gauged 
away. After this transformation, the Hamiltonian reduces 
to % = —ihvoT^cr ■ V. Because of the chiral anomaly®^, 
however, the aforementioned gauge transformation gen¬ 
erates an additional term in the Lagrangian Cem that 
describes the coupling between light and 3D WSMs®®: 

C,^ = ^{E^-B^)-pcf + J-A + Ce, (3) 

OTT 

where 

C 0 = -^9{r,t)E-B . (4) 


Here, a = e^/{hc) ~ 1/137 is the usual QED fine- 
structure constant and 0 {r,t) = 2 {b ■ r — hot) is the so- 
called axion angle. The additional axion term Cg changes 
two of the four Maxwell equations, i.e.®® 

V-£; = 4^(p+^6.b) (5) 


and 


1 dE 
c dt 


-b V X B = 


c V 


a 

27r^ 


cb X E 



( 6 ) 

Earaday’s law, x E = —c ^dBjdt, and the equation 
stating the absence of free magnetic poles, V • B = 0, are 
unchanged. 
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Eliminating B, we obtain the following modified wave 
equation in a WSM: 


1 d^E _ 47r dJ 

c2^ “ X ( X ) - 

2a, dE 2a, „ ^ 

- b X —-6oV X E 

TTC at TTC 


( 7 ) 


As usual, we now need an expression that relates the 
current J to the electric field E, which we proceed to 
derive by utilizing a semiclassical approach. 
Semiclassical Boltzmann transport theory in WSMs. —In 
the linear response regime®^®, the Fourier components 
J{q,uj) of the induced current density are linearly de¬ 
pendent on the total electric field (i.e. the sum of the ex¬ 
ternal field and the Hartree contribution), i.e. w) = 

where is the optical 

conductivity tensor and the indices to, n run over the 
Cartesian coordinates x, y, and 2 . We will work in 
the local approximation®, which is justified in the limit 
q where Rc = ud/wc is the cyclotron radius, with 

ujc = eB/{mcc) the cyclotron frequency, TOc = £f/v‘^ 
the WSM cyclotron mass, and ep the Fermi energy. We 
therefore have arnn{q,xi) « a-mn( 0 ,a;) = ^mnCTb(w) + 
o’mra(w), where crb(w) = —i(eb — l)a;/( 47 r) is the bound- 
charge contribution, while (Tmn(w) represents the free- 
charge contribution^®. The latter quantity can be cal¬ 
culated by utilizing a semiclassical Boltzmann transport 
approach, which is justified when uj,u!c £f/S. Also, 
we focus on the collisionless ujt ^ 1 regime, where 
T = min(rintra, Winter) Is the shortest between intra-node 
biitra (s-g- due momentum-non-conserving collisions) and 
inter-node Ti„ter scattering times. 

Our interest in this work in on doped Weyl semimetals 
(/cbT <§; ep, where T is temperature). The situation 
for k^T ^ £p (neutral Weyl semimetals) is much more 
complicated as one needs to include finite temperature 
effects (thermally excited carriers) and disorder. 

For a given chirality y = ± of a single Weyl node, the 
semiclassical Boltzmann equation (SBE) reads as follow- 

y24,25. 


mg- 


df 

^+P 9 -Vp/g + rg-V./g = 0 . 


( 8 ) 


Here, fg is the electron distribution function. In the pres¬ 
ence of a static magnetic field B and a time varying elec¬ 
tric field E, the semiclassical equations of motion are^® 


r = Vg{p) - p X ng{p) (9) 

and 

p = —^E — X B . (10) 

h Tic 

The first term on the right-hand-side of Eq. (9) is Vg{jp) = 
h~^\7peg{p), defined in terms of an effective band dis¬ 
persion £g{p). In topological metals such as WSMs, 
this quantity acquires a term due to the intrinsic or¬ 
bital moment^®, i.e. £g{p) = £o{p) — 'fT^g{p) ' where 
£o{p) = RvdP with p = IpI is the ordinary conduc¬ 
tion band energy while mg{p) is the orbital moment^®, 
i.e. mg{p) = 'yeeo{p)Ctg{p)/{hc). Here, 7 is a dimen¬ 
sionless control parameter and ^g{p) = —gp/{2p^) is 


the WSM Berry curvature^®. The parameter 7 takes two 
values: 7 = 1 is what one should use, while 7 = 0 is 
what one should use to artificially discard the impact of 
the orbital magnetic moment. 

Using Eqs. (9)-(10) and carrying out straightforward 
algebraic manipulations we find: 


r = rB + rE , 


( 11 ) 


where 


rB = Dg (P) ^Vg{p) + —ng{p) ■ Vg{p)B 


and 


t-e = Dg\p)^E X ilgip) . 


( 12 ) 


(13) 


Here, Dg{p) = [1 -I- eVLg ■ B/{hc)\ and the group velocity 
Vg{p) is given by 


Vg{p) = vuP 




( 14 ) 


with p — p/p and Hg(p) = |rig(p)|. Similarly, we find: 
P = Pb + Pe with pb = -D-^{p)eVg{p) x B/c and 
Pe = -D-^{p) [eE + e{E ■ B)ng{p)/{hc)]. 

Let us start by setting E — 0 while keeping B = Bz 
finite. In this case, the SBE ( 8 ) is solved by 


f^°Hp) 


exp 


£g(p) -gp 

kBT 


(15) 


when in the collision integral Ig we take /eq = (p ). 

We now want to solve the SBE up to first order in the am¬ 
plitude of a homogeneous time-dependent electric field, 
E = .E(w)e“*“*. To this end, it is useful^^ to exploit 
the symmetry of system by using cylindrical coordinates: 
p = {yjp^ — pI cos{ip), — pI sin(:/ 3 ),p^)®". We seek a 
solution of the SBE of the form 


fg{p,t) = f^°Hp) + Sfg{p,t) , (16) 


where Sfg{p,t) is linear in E and is parametrized as fol¬ 
lowing: 


Sfg{p,t) 


djf^ 

deg 


(A_e*^ -k A+e-*‘g -k Aq) , 

(17) 


with X±^o = i^±,oiPTPz)- Tbe linearization of the SBE 
( 8 ) is greatly simplified by the observation that p^ and p 
are constants of the motion in the limit E ^ 0. 
Inserting (I6)-(I7) in Eq. ( 8 ) we find 


X± = evuS 


^^hcp^ Vp^^,E,±tEy 


1 ^ u 

^2hcp^^ 


2p i{uj±ujf) 


(18) 


and 
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-^0 = ei'D 






( ® 

,2 H2' 

\2hc2 

pi 


1 ^ T3 

^2hcp3^ 


Pz 


ILO 


(19) 


with 

1 ^ P^ u 

1 - 19t — 

^l = ^t{p.Pz)=^c ■ ( 20 ) 

In writing Eqs. (18)-(19) we have introduced another di¬ 
mensionless control parameter, 6. This takes two values: 
i5 = 1 is what one should use, while 5 = 0 is what one 
should use to artificially discard the anisotropy in the 
distribution function. 

The distribution function determines the total current 
carried by electrons at each Weyl node: 

^9 = -^ J-^Bgip)rfg . ( 21 ) 


The factor Dg{p) ensures that the number of states in the 
volume element remains constant in time^®. For the sake 
of convenience, we decompose the current density per 
node in the sum of three terms: Jg = -b -I- 6Jg, 
where Jg^l = —e(27r)“^ / d^pDg{p)r\fg^^ with X = E,B 

and 6Jg = -e(27r)-^ f (PpDg(p)rB6fg. 

We first examine in detail the dependence of Jg^l on 
the static magnetic field B. In the weak magnetic field 
limit, we expand the distribution function fg^'^ in powers 
of B, up to second order: 


o py 

f!,%9iP)] « /f [^o(p)] - ^ 


(0) 


l-rrigip) ■ B 


-b - 


1 52 /i°’ 


2 dej 


£o(p) 
[7mg(p) • B ]2 . 


£o(p) 


Taking the limit T <C Ef/^b we find 


( 22 ) 


j{Q) _ e^pp 

g,E ^247r2fi,ceF 


B X £:(a;)e-*‘^‘ 


(23) 


and 


— Q 

-Jg.B - 9 


e^e-p 

Att'^ItPc 


B 


(24) 


We see that Jg% is a) independent of the chirality g of 
the Weyl node, and therefore leads to a finite correction 
to the ordinary Hall conductivity—see Eq. (29) below— 
and b) proportional to the dimensionless parameter 7 . 
Because of b), Eq. (23) originates from the orbital mag¬ 
netic moment mg(p). The term in Eq. (24) is pro¬ 
portional to the chirality g of the Weyl node, and there¬ 
fore has no effect on the total current J = Ylg Jg but 


yields a finite axial current Jax = J^gPJg^ agreement 
with Refs. 28 and 29. 

We then evaluate the quantity 6Jg and obtain the op¬ 
tical conductivity tensor We first consider SJg^^- 

By retaining all terms of second order in the ratio Swc/Mj 
we find^° 

V 

crzz{uj)=i —, (25) 

TTUJ 

where V = ’Ke^n^/mc is the Drude weight, rig = 
ep/{3Tr‘^h^v^) is the electron density, and 

^ ^ , 3 ;i 2 2 ^ 3 ^ 7[27(5 + 5) + 1W-25] 

= ^ + + ^ -15-— 

(26) 

Eqs. (25)-(26) are the most important results of this Sec¬ 
tion. Setting i5 = 7 = 1, we obtain the desired result for 
the longitudinal conductivity in the presence of a weak 
magnetic field: 




(27) 


Because of the non-trivial dependence of lu* on p in 
Eq. (20), the calculation of SJg^x, ^Jg,y at an arbitrary 
frequency w is not straightforward. This calculation, 
however, notably simplifies in the low-frequency w <C Wc 
limit, which is relevant for helicons. In this limit and 
after setting 7 = 5 = 1 , we find 


and 


1 

20 £2 j > 


rwc y 20 j 


, Voj 


'KUJf. 


1 - 


(28) 


(29) 


The remaining off-diagonal elements of the optical con¬ 
ductivity tensor, such as o’yz, etc., vanish identically 
for symmetry reasons, independently of the frequency w. 
Helicons and plasmons in WSMs. —Using the wave equa¬ 
tion (7) and the semiclassical result for the optical con¬ 
ductivity tensor crmn(w), we seek for collective modes of 
doped WSMs subject to a weak static magnetic field. To 
this aim, it is useful to introduce the dielectric tensor 


— ^£m^h “b 


47rf 

UJ 


^ Cm ^£mn 


ac 

27r2 



(30) 


and Mem = c?{q'^5em - didm) - In Eq. (30), eemn 

is the 3D completely antisymmetric tensor and the Latin 
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indices £, m, and n run over the Cartesian coordinates x, 
y, and z. Finally, a sum over n is intended. 

The zeroes of the determinant of A4 correspond to the 
self-sustained modes of a doped WSM. Following stan¬ 
dard practice®, we focus on two special cases: i) q par¬ 
allel to the static magnetic field B, i.e. q — qz and ii) 
q orthogonal to B. When q = qz, four collective modes 
appear: three gapped modes, which are characterized by 
an energy of the order of Fermi energy, and a gapless 
mode, the helicon. If the wave vector q is orthogonal 
to B, we find only the three gapped modes, while the 
helicon solution is absent. 

After straightforward algebraic manipulations, we find 
the helicon dispersion relation in the long-wavelength 
limit: 


^h(9 0) 


2abocq/'K + c^q^ 
uj'^/uJc + 2acbz/TT 


( 31 ) 


where Wp = dTrUee^/rUc is the 3D plasma frequency in a 
WSM. Eq. (31) is the most important result of this Arti¬ 
cle. Note that Eq. (31) is independent of the background 
dielectric constant Due to the time-like component 
bo of the axion angle, the helicon dispersion relation in 
a WSM is linear in q rather than quadratic, the latter 
functional dependence on q being the one occurring in 
ordinary metals—see Eq. (1). Even for 5o = 0, the he¬ 
licon frequency differs from the textbook result (1), in 
that the effective helicon mass mh in a WSM depends 
on the component of b along the direction of the static 
magnetic field B: ruh = (2 lOc(?') -f ahbz/{TTc). 

Before concluding, we comment on the gapped collec¬ 
tive modes. For the sake of simplicity, we set B = 0 in 
the following analysis. In the long-wavelength limit, we 
find that the three gapped modes np^\{q) with A = 1, 2, 3 
are given by: 


Dp.i(g = 0)=w_ 

Dp.2(<7 = 0) = , (32) 

Dp.3(g = 0) = a;+ 


where uj± = yj(acb)'^/(ireh)^ -l-Wp/eb ± acb/{Tr€h), with 

b = |6|. Very interestingly, we find that, unlike in an 
ordinary non-topological metal®, the degeneracy of the 
three gapped collective modes at g = 0 is lifted by the 
presence of the axion term Cg in the electromagnetic re¬ 
sponse. This is due to the fact that WSMs are opti¬ 
cally gyrotropic media®® with gyrotropy parameter pro¬ 
portional to b = |6|. Since the energy of the gapped 
collective modes is comparable to the Fermi energy, an 
accurate description of these modes requires the inclu¬ 
sion of the inter-band contribution (T^^‘®''(a;)®'^^®® to the 
optical response, which has been neglected so far in our 
semiclassical approach: 


winter/, 

<^u (^) 


acLO 

127rwD 


<d{huj 


2eF) 



4A2 
4e| — 

( 33 ) 
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FIG. 1. Dispersion relations hi}p^x{q) (in meV) of collec¬ 
tive gapped (“plasmon”) modes in doped Weyl semimetals as 
functions of wave vector q, in units of kp = ep/ilivu)- Such 
dispersions have been obtained by i) generalizing Eq. (32) to 
finite q kp and ii) including inter-band effects, as from 
Eq. (33). Dashed line: A = 1. Solid line: A = 2. Dash-dotted 
line: A = 3. Panel (a): results for q parallel to b. Panel (b); 
results for q orthogonal to b. Results in this figure have been 
obtained by setting £b = 5, £f = 40 meV, Vp, = c/1000, 
bo = 0, b = (0.01,0.01,0.01)7r/a, and A = kvor'/a with 
a = 3.5 A. This choice of microscopic parameters is justi¬ 
fied by recent experimental results in NbAs and TaAs^’^^^. 


Here, I = x,y,z and A is an ultraviolet cut-off. Fig. 1 
shows the dispersion relations r2p_A(?) of the A = 
1,2,3 gapped collective modes, as calculated by adding 
cr^^*®''(a;) to the intra-band semiclassical contribution 
(^xx(!^) = o'yyi'-^) = crzz(w), i.e. Eq. (27) evaluated at 
B — 0 and for j = 6 = 1. In Eig. 1(a) [(b)] the wave vec¬ 
tor q is parallel [orthogonal] to the space-like component 
b of the axion angle. The main effect of the inter-band 
contribution on the collective modes is to redshift their 
gaps at g = 0. This can be easily explained by recog¬ 
nizing that, in the long-wavelength limit, the inter-band 
contribution (w) to the optical response can be de¬ 
scribed, to a very good approximation, as a renormal¬ 
ization of the background dielectric constant, i.e. Cb —t 
eb + aclog[|4A2/(4£| - S2u;2/eb)|]/(37ruD). In Eig. 1(a) 
we find two transverse modes with a quadratic disper¬ 
sion relation and a dispersionless longitudinal mode. We 
note that in this case there is no mixing of the two trans¬ 
verse modes with the longitudinal mode, exactly like in 
an ordinary metal. What is peculariar to WSMs is that 
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when q is tilted away from 6, the lowest-energy trans¬ 
verse mode hybridizes with the longitudinal mode. This 
effect is maximal for q orthogonal to b, as in Fig. 1(b). 

In summary, we have evaluated the optical conduc¬ 
tivity tensor of a 3D Weyl semimetal from semiclassical 
Boltzmann transport theory, with the inclusion of the 
orbital moment mg{p) and anisotropic contributions to 
the distribution function. A general expression for the 
longitudinal conductivity is reported in Eq. (25). We 
have used the calculated optical conductivity tensor to¬ 
gether with the axion contribution (4) to the standard 
electromagnetic Lagrangian to find the collective modes 
of a 3D Weyl semimetal. We have demonstrated that 
the axion term dramatically alters the helicon dispersion, 


Eq. (31), with respect to that in non-topological metals, 
Eq. (1). Finally, we have highlighted axion anomalies 
in the gapped sector of collective excitations. Fig. 1, by 
taking into account inter-band corrections to the semi¬ 
classical (intra-band) optical response. 
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In the spirit of Drude theory, the dissipative part of the 
conductivity tensor can be obtained by replacing oj —>■ 
oj + i/r in Eq. (25), where r is a phenomenological re¬ 
laxation time. We find (w) = CzzCrnjiA — iair), where 
(td = ^n^T/m^ is the usual Drude formula. We are now in 
the position to make contact with the existing dc trans¬ 
port literature®"^’®^. We first notice that the authors of 
Refs. 24 and 31 neglected corrections due to the intrinsic 
orbital moment mg{p) and terms that are responsible for 
anisotropies in the distribution function: we therefore set 
(5 = 7 = 0 in Eq. (26). We then take the dc lot —>■ 0 limit in 
the expression for cri^^cj): in this limit we recover the re¬ 
sult of Refs. 24 and 31, i.e. aiV (0) = e*v^B^T/{4:'R^hc^ep). 
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We note that even for 7 = 5 = 1, the quantity (0) re¬ 
mains an increasing function of , as in Refs. 14 and 19. 
Finally, as a caveat, we mention that one should exercise 
some care in using Drude theory in topological materials, 
since the latter fails to capture side jump and skew scat¬ 
tering contributions to transport: see, e.g.. Ref. 32. 
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